Explanation and observability of diffraction in time 
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Diffraction in time (DIT) is a fundamental phenomenon in quantum dynamics due to time- 
dependent obstacles and slits. It is formally analogous to diffraction of light, and is expected to play 
an increasing role to design coherent matter wave sources, as in the atom laser, to analyze time- 
of- flight information and emission from ultrafast pulsed excitations, and in applications of coherent 
matter waves in integrated atom-optical circuits. We demonstrate that DIT emerges robustly in 
quantum waves emitted by an exponentially decaying source and provide a simple explanation of 
the phenomenon, as an interference of two characteristic velocities. This allows for its controllability 
and optimization. 
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Diffraction in time is a fundamental quantum dynami- 
cal effect first studied by Moshinsky [1]. ID matter waves 
released through a time-modulated aperture or encoun- 
tering a time-dependent obstacle (for 2D and 3D cases 
see 0, El ) show temporal quantum penumbras and inter- 
ference patterns similar to the diffraction of light behind 
spatial slits and obstacles. Understanding and control- 
ling DIT is more and more relevant as a result of the in- 
creasing manipulability of coherent matter waves, in par- 
ticular in ultracold atomic gases and/or with ultrashort 
laser pulses. DIT will affect, for example, the intended 
applications of atom lasers [4| , dynamics of matter waves 
emitted by ultrashort laser excitations [H[ , matter- wave 
circuits [6[, and time-of- flight techniques [7[. DIT may 
also lead to temporal versions of diffractometers, grating 
spectrometry, and holography. 

The original and most studied setting for DIT is the 
Moshinsky shutter (MS). It consists in a sudden release, 
by opening a shutter, of a semi-infinite plane- wave beam 
characterized by a "carrier" velocity. The particle den- 
sity as a function of time at an observation point is for- 
mally analogous to spatial Fresnel diffraction by a sharp 
edge If the shutter, when closed, has reflecting am- 
plitude R = 1, the same results are obtained from a point 
source with a sharp onset and constant emission there- 
after [§]. Many works have applied and modified MS to 
study different quantum transients and, adding a poten- 
tial, resonance scattering, buildup and decay, and tun- 
neling dynamics, see e.g. [{J |To| and reviews in d, 
Experimentally, a DIT oscillatory pattern was first ob- 
served by Dalibard and coworkers with cold atoms falling 
by gravity and bouncing off a mirror made of evanescent 
light that could be switched on and off [12]. DIT through 
a related time-energy relation has been observed for cold 
neutron experiments too [7[ . Interferences from two time 
slits and time-analogues of diffraction from a grating have 
been described for cold atoms [H, [Hj and ionizing atoms 
with ultrashort laser pulses [gj . There are also analogs of 
the original MS in the field of coherent transients due to 



frequency-chirped weak lasers [HI . 

DIT may be suppressed or averaged out by apodiza- 
tion, noise and decoherence, or unsharp carrier velocity 
distributions [HEl, so tne observability of MS-DIT with 
matter waves has been considered a difficult task [8[ . We 
shall see, however, that the effect is rather robust and oc- 
curs quite generally in waves emitted by an exponentially 
decaying resonance. 

A second problematic aspect of MS-DIT is the lack of 
a simple and intuitive understanding of the phenomenon. 
The usual, geometrical "explanation" in terms of a Cornu 
spiral []], [2, 8] does not provide a simple physical picture 
although some insight is gained by its construction via 
Fresnel time-zones and Huygens principle, as in spatial 
diffraction [2[. Also, an attempt was made in [16| to seek 
an explanation in terms of the Wigner distribution but, 
as recognized by the authors, the interpretation of the 
results remained ambiguous due to the lack of positivity 
of the Wigner function. 

In this paper we address the observability and interpre- 
tation of DIT, as well as some important consequences. 
They are linked to each other since a simple physical ex- 
planation of DIT will also provide the key to observe and 
control it. The starting point is the realization that sys- 
tems that decay exponentially due to a resonance, such 
as cold atoms in magnetic or optical traps [ItJ escaping 
from their initial confinement, may show DIT at a dis- 
tance from the trap. The density or flux oscillations will 
be identified as an interference effect characterized quan- 
titatively with a simple analytical model [18|. We shall 
thus be able to predict and design optimal conditions for 
its observability, and treat on the same footing the stan- 
dard constant emission after a sharp onset, and the expo- 
nentially decaying source, by modifying continuously the 
imaginary part of the emission pole. Figure [TJ discussed 
later in more detail, shows the unnormalized density at 
an observation point away from the source. The upper 
curve corresponds to ordinary MS-DIT oscillations. The 
amplitude of the oscillations at the observation point de- 
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FIG. 1: (Color online) Unnormalized density versus time at 
x = 1000 for a constant or exponentially decaying source. 



creases with time, and their frequency depends on time, 
tending to a constant. The other two curves correspond 
to exponentially decaying sources with different lifetimes. 
The oscillations are essentially the same as in the stan- 
dard MS, modulated by the exponential decay. 

The exponentially decaying source model. We shall use 
a model that captures the essence of resonance decay 
from a trap and describes analytically the external wave 
function without the complications and peculiarities of 
particular confinements [18j. We adopt the same nota- 
tion as in [18] with dimensionless position x, time t, and 
wave function ip obeying formally a Schrodinger equa- 
tion for a particle of mass 1/2 and h= 1, idip(x,t)/dt = 
—d 2 ifj(x,t)/dx 2 . The unit of length corresponds to the 
inverse of the real part of the carrier wavenumber and 
the unit of time to the carrier period divided by 2tt. The 
complex dimensionless wavenumber fco = fcoi? + ^o/ an d 
frequency of the carrier ujo = ujor + zc^oj, obey the dis- 
persion relation ujq = fc 2 , = (1 + ifco/) 2 , so ujor = 1 — fc 2 ^ 
and uoi = 2fcoi, with koi < and koR = 1. The exact 
unnormalized solution to the Schrodinger equation for 
the free particle subjected to the source boundary con- 



dition i/j(0,t) = e~ iuJot Q(t), ujor > 0, uj oi < 0, can be 
constructed by a superposition of plane waves. The re- 
sulting integral can be expressed in terms of known func- 



tions, ip(x,t) 
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where 



w(z) := e~ z erfc(-iz), ^ = ±(1 + i)y/t/2k (l =F r/t), 
and k s = x/2t, r — x/2ko are a "saddle point" wave 
number and a complex characteristic time. For an ob- 
servation point x, the saddle velocity is time dependent, 
v s = 2k s = x/t. Figure [U shows the unnormalized den- 
sity \i(j(x,t)\ 2 for different koi to illustrate the essential 
continuity of oscillation phenomena when varying fco/. If 
one particle is emitted, the normalized wave function is 
$(x,t) = [f™dt J(0,t)] _1 / 2 ^(a:,t), where J(x,t) is the 
dimensionless flux. 

The essence of BIT. The wavefunction -0, for times 
shorter and larger than |r| [18], can be accurately 
approximated by contributions from the two criti- 
cal points of its defining integral, saddle and pole, 



FIG. 2: (Color online) Normalized probability densities versus 
time for fco = 1 — 0.0015i at x = 1000. Exact (red solid line), 
approximate (Eq. (pQ), black dashed line), saddle term (blue 
circles), pole term (orange triangles), and interference term 
(Eq. ([2]), green squares). 
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As a result of contour deformation along the steepest 
descent path from the saddle, the pole term contributes 
from the time when Im(^Q + ^) = 0, t c = x/[2(l+fc /)], and 
decays exponentially thereafter. In a pictorial, classical 
association [19] , the particle arriving at (x,t) with veloc- 
ity vq = 2 must have been released at a time x/vq from 
the source which emits particles exponentially. The sad- 
dle velocity x/t is the one required for a classical particle 
released from (0,0) to arrive at (x,t). Saddle trajecto- 
ries may thus be pictured as the result of a burst or "big 
bang" emerging from the source with all possible veloc- 
ities at t = 0. These classical pictures are useful but, 
unlike long-time deviations from exponential decay [l9| . 
DIT cannot be explained with them alone. It is a quan- 
tum interference phenomenon as shown by the structure 
of the unnormalized density 



\^{x,t)\ 2 = \Mx,t)\ 2 + \Mx,t)\ 2 e[lm(u+)} 
+ 2Re[Mx,tWo(x,t)]e[Im(u+)], 

where the asterisk denotes complex conjugation, 
interference term is 
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Light does not show DIT in vacuum because 
there is no dispersion and no interference of this kind. 

Figure [2] shows the agreement at times larger and 
shorter than |r| between exact and approximate wave 
functions. The pole and saddle terms separately do not 
oscillate in time, whereas the interference term, Eq. (j2j), 
reproduces accurately the characteristic oscillations of 
DIT. 

Characterization and observability of DIT. The fre- 
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FIG. 3: (Color online) Time intervals T n+ i — T n between 
two consecutive maxima: exact (circles), and approximation 
from Eq. J3} (squares). Same parameters as in Fig. [2 The 
symbols size is to help the eye, not related to errors. In the 
inset: nonlinear position of the first maximum versus time 
for koi = —0.08 (solid line) and onset of the pole term x c = 
2(1 + koi)t (dotted line). 



quency of the DIT oscillation depends on the interference 
of the saddle and pole frequencies k% and uor and, as k s 
depends on time, the DIT oscillation period is not con- 
stant. From Eq. (|2j) we can infer the position of the n-th 
maximum. For |fco/| <C 1, the s'mcj) term of Eq. (|2]) tends 
to vanish so the DIT oscillations are essentially due to the 
cos(/> term. The maxima correspond to (f)(x,T n ) = 2nir 
at times 



FIG. 4: (Color online) Saddle and pole terms for koi = —0.03 
(solid blue line) and koi — —0.13 (long-dashed red line), x = 
80. The inset shows the corresponding exact densities. The 
shorter lifetime suppresses DIT. 
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where n = 0,1,2, ... (n = is for the principal max- 
imum). The interval T n+ i )U = T n+i — T n between two 
consecutive maxima is in good agreement with the exact, 
numerically calculated period, see Fig. [3l The small dis- 
crepancy at n = can be attributed to the dependence 
on time of the factors multiplying cos <j) and the proximity 
of|r|. 

For large times the period of the DIT oscillations tends 
to the carrier period, lim T n+ i ?n = 2tt/ujor. The ampli- 

tude of the oscillations decays relatively slowly compared 
to the pole term, as e ^o/^-3/2^ gee q ^ut expo- 
nentially faster than the saddle term. 

According to Eq. (J3j), Tq is not a linear function of 
x. For example, in the limit koi —> the motion of the 
first maximum is described by xq = 2T — ^/3Tqtt. Thus, 
even though an asymptotic velocity may be denned, 2(1 + 
koi) in the general case, see the inset of Fig. [3j there 
is no oblique asymptote for this function. Therefore a 
naive linear extrapolation back to the origin at some large 
distance fails to provide the instant of the source onset. 
In other words, the times in which the tangents to xo(t) 
cut xo = have no definite limit, in spite of the well 
defined asymptotic velocity. This is an example of the 
importance of DIT to correct simple classical-dynamical 
extrapolation from asymptotic wave features to extract 
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FIG. 5: (Color online) Density \ij)(x,t)\ 2 for koi 



-0.003 



showing the transition from pure exponential decay to a DIT 
pattern. 



emission characteristics, as practiced e.g. in the analysis 
of ionization by ultrashort laser pulses [20| . 

In our dimensionless description two factors affect the 
visibility of the DIT pattern: the observation position 
x and the lifetime. Figure |4] shows the moduli of the 
logarithm of the pole and saddle densities for two dif- 
ferent lifetimes. The pole term is a semi-infinite straight 
line which begins when the pole is crossed by the steepest 
descent path passing along the saddle in the complex mo- 
mentum plane, at t c \ the saddle term shows a maximum 
near |r| and decays from there slowly. There may be 
up to two intersections of the two terms, one near the ar- 
rival of the main front, and one at a long time that marks 
the transition to post-exponential decay [18|. When the 
saddle and pole terms are similar or close enough the in- 
terference oscillations appear. The interference region of 
our interest here is the one following the main front be- 
cause it relates by continuity to ordinary MS-DIT in the 
limit koi — ^ 0; it is also much more easily observable than 
the oscillations at large times because of the magnitude 
of the amplitudes. 

The oscillations are evidently not present at the source 
x = 0, and will be small at small distances, x < 1, be- 
cause of the rapid decay and separation from the pole 
term of the saddle term in these conditions, see Fig. [5] 
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FIG. 6: (Color online) Stability of DIT for decay from US(x) 
(Winter's model). The initial states are the ground states of 
the two wells in the inset. At t = the right wall, at x — 0, is 
substituted by the delta. L = 3.14, x = 157.05, U = 161.35, 
V = 202.72. 

The saddle term beyond the main front arrival increases 
with x [18|- In the opposite extreme of very large x, it 
eventually dominates entirely and stays above the pole 
term at all times, suppressing DIT and even exponential 
decay [HI . Between these two extreme scenarios there is 
an ample range of x for which DIT is prominent. The 
slope of the pole term also plays a role. For larger values 
(smaller lifetimes), pole and saddle contributions sepa- 
rate more rapidly leading to fewer visible DIT oscillations 
which may actually disappear for small enough lifetimes. 

To estimate the domain where some oscillations are 
seen before the decay is too strong we may solve T^o < 
Ntq for a small TV, where To = l/4|fco/| is the lifetime. 
This gives an explicit but lengthy expression. For N = 5 
and in the koi —> limit, x < 30tq . 

From the previous discussion it might seem that a very 
long lifetime is always preferable to attain DIT. Neverthe- 
less long lifetimes also imply a weaker signal because of 
the normalization. The consequence of opposite tenden- 
cies is an optimal lifetime-position point. A good mea- 
sure of the visibility of DIT is the difference A between 
the second maximum and the previous minimum of the 
normalized probability density, see Fig. [2j The optimal 
parameters are found to be koj = —0.03, x = 60. 

Model independence of the results. We have described 
the close connection between DIT and resonance decay. 
DIT will occur when contributions from different reso- 
nances are well separated, which generally requires nar- 
row and/or strong confinement. DIT does not depend on 
the specific properties of the model used so far. We have 
checked the robustness of DIT from exponential decay 
explicitly with several additional models. Winter's decay 
model [21[ describes the decay of the ground state of the 
square well between — L and when the right infinite wall 
is substituted by a US(x) potential. The wave function 
outside the trap tends to the source model wavefunction 
for large U [H[ . Moreover DIT does not depend dramati- 



cally on the strict confinement of the initial wave function 
on a finite domain. To show this we have calculated the 
decay of the ground state of a well with a finite right wall 
once this wall is substituted by the delta. This produces 
a different fast forerunner at x, but the part associated 
with the dominant, lowest energy resonance remains es- 
sentially stable showing DIT as for the infinite wall, see 
Fig. [6l Moreover we have observed the same stability for 
finite- width barriers. DIT also survives a smooth source 
onset [4], and again may be observed after the passage 
of some onset-dependent transients. As for the effect of 
collisions, in the mean field regime DIT is enhanced for 
attractive interactions [8|- 

Let us finally point out the possibility to observe DIT 
in periodic structures [9] such as optical lattices, or other 
physical systems that realize a tight-binding model, for 
example periodic waveguide arrays that provide a classi- 
cal, electric field analog of a quantum system with expo- 
nential decay |23| . 
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